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In this paper we give an exact infinite-series expression for the bi-partite entanglement entropy 
of the quantum Ising model both with a boundary magnetic field and in infinite volume. This 
generalizes and extends previous results involving the present authors for the bi-partite entan- 
glement entropy of integrable quantum field theories, which exploited the generalization of the 
form factor program to branch-point twist fields. In the boundary case, we isolate in a universal 
way the part of the entanglement entropy which is related to the boundary entropy introduced 

{Sj . by Affleck and Ludwig, and explain how this relation should hold in more general QFT models. 

We provide several consistency checks for the validity of our form factor results, notably, the 
identification of the leading ultraviolet behaviour both of the entanglement entropy and of the 

"V* ■ two-point function of twist fields in the bulk theory, to a great degree of precision by including 

up to 500 form factor contributions. 
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1 Introduction 



Entanglement is a fundamental characteristic of quantum systems, which has great importance, 
for instance, in the context of quantum computing. Any measure of entanglement is likely to give 
a good description of the quantum nature of a ground state. Many measures of entanglement 
have been devised, see e.g. PQ-[S]. However, perhaps due to its geometric character which makes 
it more theoretically tractable in systems with many degrees of freedom, the entanglement 
entropy [T] has attracted great interest in theoretical physics. The entanglement entropy can 
also be argued to be a better measure of fundamental properties of the ground state than 
correlation functions, as it is not associated to a particular observable, but rather to a sector of 
mutually local observables. In this paper we continue the research initiated in [BJ [7J |8] on the 
bi-partite entanglement entropy in one-dimensional quantum models with many local degrees of 
freedom, focusing on the effect of boundaries. 

The entanglement entropy is a measure of quantum entanglement between the degrees of 
freedom of two regions, A and its complement A, in some quantum state Consider a quantum 
system, with Hilbert space TL = TLa <S> 7~(-a-> in a pure state The bipartite entanglement 
entropy Sa is the von Neumann entropy associated to the reduced density matrix p A of the 
subsystem A, 

p A = Trn A m^\), (1.1) 

given by 

S A = -Tr HA (p A logp A ) = - lim ^-Ti HA (p n A ). (1.2) 

n-*i an 

The expression with the n-limit on the right-hand side is often referred to as the replica-trick. 

Let us consider now the scaling limit of the quantum system, describing the universal part 
of the behaviour near a quantum critical point. It is obtained by approaching the critical point 
while letting the length of the region A go to infinity in a fixed proportion with the correlation 
lengths. The result is a quantum field theory (QFT) model, which we will take throughout 
to possess (1 + 1-dimensional) Poincare invariance. In what follows, we consider only the case 
where A is a connected region. 

In the QFT context, the expression on the right-hand side of the second equation of (|1.2|) 
can be understood, for n a natural number, as a normalised partition function for the model 
on a Riemann surface with two branch points, with n sheets cyclicly connected, or on a surface 
with two conical singularities of angles 2im (9j fTTJl [TT] (see also the explanations in [SI [7j ) . There 
is only one way of associating such branch points to well-defined local QFT fields. This was first 
done in [6] in the case without boundaries (with A a region in the bulk). There, it was shown how 
to relate the entanglement entropy in two-dimensional QFT with a two-point function of certain 
local fields defined in a model consisting of n copies of the original model, called branch-point 
twist fields (section [2]) : 

S h / lk (r) = - lim ^-Z n e 4A "(0\f( Xl )T(x 2 )\0}. (1.3) 
n-+l dn 

Here A has length r = \x\ — x%\ and (0| • • • |0) denote correlation functions in the n-copy model; 
the state |0) is the vacuum state of the latter. The derivative with respect to n involves an 
appropriate analytic continuation in n of the correlation function, which is assumed to be in 
correspondence with the conical-singularity interpretation. We will not discuss further in the 
present paper the subtleties and assumptions involved in this analytic continuation - see the 
discussion in [7] for more details. The constant Z n , with Z\ = 1, is an ra-dependent non-universal 
constant, e is a short-distance cut-off which is chosen so that dZ n /dn = and, finally, A n is the 
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conformal dimension of the counter parts of the fields T, T in the underlying n-copy conformal 
field theory, 

A„ = ^(„-i), (1.4, 

which can be obtained by CFT arguments [111 |6j and where c is the central charge. 

In [HIE!, the two-point function (0\T(xi)T(x2)\0) was studied at large distances r for all 
1+1-dimensional integrable QFTs on the line. The main feature of these models is that there 
is no particle production in any scattering process and the scattering {S) matrix factorizes into 
products of 2-particle S-matrices which can be calculated exactly (for reviews see e.g. [12]- [16]). 
Taking the S-matrix as input it is then possible to compute the matrix elements of local operators 
(also called form factors) . This is done by solving a set of consistency equations |17l [T8] . This 
is known as the form factor bootstrap program for integrable QFTs. In [6j[7], this program was 
used and generalised in order to compute the two-particle approximation of (0|T(xi)T(x2)|0), 
and to obtain the leading correction to saturation of the entanglement entropy. This leading 
correction was observed to be very universal, as it is independent of the scattering matrix. 
In fact, by similar techniques, this universal correction was also observed to hold outside of 
integrability [5]. 

In this paper we generalise the construction above in order to include the presence of one 
integrable boundary. The study of integrable QFTs with boundaries has attracted a lot of 
attention in the last two decades (see e.g. |19|-|24j). The present work, will make extensive use 
of the results of S. Ghoshal and A. B. Zamolodchikov [22], particularly the explicit realization of 
the boundary which they proposed. Their work provided also a detailed study of the Ising model, 
for which the integrable boundary conditions were classified and the corresponding reflection 
amplitudes computed. These reflection amplitudes will provide a crucial input for our entropy 
computation. 

Let us therefore consider a family of integrable boundaries parametrised by the dimensionless 
constant 

k = 1 - h 2 /(2m) e (-co, 1) , (1.5) 

related to a uniform magnetic field h affecting the boundary Ising spins. We study the entangle- 
ment entropy between a region A that extends from the boundary to a distance r, and the rest. 
For k < 0, we obtain the full large-distance series expansion; the result depends on the reflection 
matrix at all orders. It turns out that some of the techniques necessary to obtain this result are 
also useful in the bulk case, so that as a by-product, we obtain the equivalent expansion in the 
Ising model without boundaries, extending our previous work [6J. We note that this extension 
to higher particle contributions involves subtleties that were not present in the two-particle case. 

We then evaluate from the form factor expansion the exact universal constant V(n) that 
relates the large-distance value of the entanglement entropy to the short-distance logarithmic 
behaviour in the boundary case, defined by 

/ -log(2r/e) + V(K) + o(l) e <C r < mT x 

6 c U , (1-6) 

— log (em) H h 0((rm) °°) r 3> m 

6 2 

where m is the mass of the particle of the Ising model and c = 1/2 is its central charge (and 
o(l) is in terms of the small combination rm). Note that the corrections to the large-distance 
saturation are exponential. The constant U is the universal saturation constant that occurs in 



5 boundary (r) 
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the bulk case, calculated in the Ising model in [6j: 




- log(r/e) + o(l) e<^r<^m 

O 

— log(em) + U + 0((rm)~°°) r > mT x 



-1 



(1.7) 



Equations ()1.7p and (|1.6f) provide universal definitions for U and V(k). Note that e here is in 
general a different short-distance cut-off in the bulk and boundary cases. The choice of these 
definitions will become clear later. 

Our main findings are 1) the observation of non-monotonicity of the entanglement entropy 
in the boundary Ising model for boundary magnetic field lower than a critical value, and 2) the 
relation between the constant V(n) and the boundary degeneracy g of Affleck and Ludwig [25]. 
Exact re-summations of form factors and CFT arguments strongly suggest the following result: 



In more general cases of massive QFT, for h associated to a relevant boundary perturbation, we 
argue that 



where s = logg is the boundary entropy in the UV (k > —00) or infrared (k = —00), and C 2 is 
the fraction of the massive ground state degeneracy that is broken by the field h. This provides 
a way to extract the boundary entropy solely from universal entanglement entropies. 

The paper is organised as follows: in section [2] we introduce the bi-partite entanglement 
entropy of a general massive boundary QFT. We establish a relationship between this quantity 
and the derivative at n = 1 of the boundary one-point function of a twist field associated to 
the QFT constructed as n non-interacting copies of the original model, similarly as discussed 
for the bulk case in (6j \7\ [8]. We provide a CFT analysis of the constant V(k) introduced 
above. We then specialize these results to the case of integrable models, introducing the form 
factor expansion for the boundary entanglement entropy. Finally we introduce the Ising model, 
for which we review the different kinds of integrable boundary conditions and the associated 
reflection matrices obtained in [22]. In section 3 we start by reviewing the form factor approach 
for branch-point twist fields and give a closed expression for all non-vanishing form factors of 
the Ising model. We proceed to a detailed analysis of the boundary entanglement entropy in 
the Ising model, giving a closed expression for all form factor contributions. The derivation of 
these formulae involves a complicated analytic continuation on the variable n, as the derivative 
is taken, followed by the n — * 1 limit. In section 4 a similar analysis is performed for the bulk 
theory, extending the results of [BJ. In section 5 we provide a analytical and numerical study 
of the ultraviolet behaviour of the entanglement entropy of the bulk and boundary Ising model, 
and evaluate V(k). In section 6 we provide a discussion of the main results, and in section 
7 we summarize our main conclusions and open problems. We provide three appendices: In 
appendix A we give a proof of several formulae which we have used for the computation of the 
entanglement entropy. In appendix B we provide alternative formulae for the individual form 
factor contributions to the bulk and boundary entanglement entropy which are more suitable 
for numerical computations. In appendix C we provide a detailed analysis of the UV behaviour 
of the two-point function of the twist fields in the bulk Ising model and extract the coefficient 
of the logarithmic term with great precision. 





V(k) = s-logC, 
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2 Entanglement entropy of two-dimensional QFTs with bound- 
aries 



2.1 General considerations 



We consider here a general massive boundary QFT with one of the particle masses being m, 
characterising the scale of all other masses, and a boundary parameter h, associated to a relevant 
perturbation of a conformal boundary condition. More precisely, we may write the action as 

S = S hulk (m) + h J dt<p(t), (2.1) 

where <ft(t) is a boundary field with dimension less than 1 and t is time. For later convenience, 
we will use k = 1 — h 2 /(2m) E (— oo, 1). 

The main idea of [6j in order to evaluate the entanglement entropy was to use n independent 
copies of the original model, which preserves integrability. The introduction of extra copies 
allows for the presence of a 7L n symmetry, associated to which the branch-point twist fields 
T(x) and T(x) can be defined. It is these twist fields that implement the non-trivial connection 
between sheets in the Riemann surface with branch points, used to calculate the entanglement 
entropy. Below we recall some of their main properties, but for more details about branch-point 
twist fields and the construction in the bulk case, see the discussions in [HIE]- 

The branch-point twist fields can be characterized as follows: let ^i, ... , \l/ n be any fields 
belonging to each copy of the original modefl- Then the equal time (x° = y°) exchange relations 
between T(x) and ^Si{y), ... , ^ n (y) can be written in the following formj]: 

%(y)T(x) = T(x)V i+1 (y) x 1 > y\ 

^(y)T(x) = r(x)*i(y) x 1 < y\ (2.2) 
for i = 1, . . . , n and where we identify the indices n + i = i and similarly 

*i(l/)T(x) = f(x)^i^(y) x 1 >y 1 , 

*i{y)f{x) = f(x)%(y) x l <y\ (2.3) 

that is T = T"'. These exchange relations indicate that a branch cut originates from branch- 
point twist fields in (Euclidean) correlation functions with insertion of fields The definition 
of the branch-point twist field is completed by saying that they are spinless, invariant under 
all symmetries of the original model, and of lowest dimension, and by specifying the CFT 
normalisation T(x)T(0) ~ |x| _4An as |x| — > + (space-like), with the conformal dimension given 
in (fOl) . Note that as n ->■ 1, T -* 1. 

We now consider the presence of a boundary. We will place the boundary at the origin of 
space x 1 = and therefore have a QFT defined on the (positive) half-line. In order to make a 
clear connection with the bulk results, we first consider the entanglement entropy S(ri,r2) in 
the case where the region A is a bulk region, extending from r\ > to ri > r%, region A being 
the rest, composed of two disconnected components, from to r% and from r2 to oo. Arguments 
entirely similar to those of [6j [7] show that the trace in (jl.2p becomes a normalised correlation 



"The vector space of fields of the multi-copy model contains the n-fold tensor product of the vector space of 
fields of the original model. Fields belonging to a copy i have the identity field 1 for all factors corresponding to 
copies j ^= i. 

^Here we employ the standard notation in Minkowski space-time: x v with v = 0, 1, with x° being the time 
coordinate and x 1 being the position coordinate. 
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function of twist fields as in (jl.3p . but with the ground state of the model on the half- line, 
\0)b- Note in particular that the boundary condition is Z„-invariant, so that the branch cuts 
originating from the twist fields can be deformed through the boundary. 
Then, the entanglement entropy is given by: 

S A (n,r 2 ) = - lim ^-Z n e 4A - B (0\f(r 2 )T( n )\0) B (2.4) 
n— >1 an 

Again, the non-universal constant Z n satisfies Z\ = 1 and dZ n /dn = 0, e is a short-distance 
cut-off, and A n is given in (II. 4h . In this equation, we have a vacuum correlation function in the 
model on the half-line, the twist-fields being at positions r\ and r 2 on the half-line. 

There are two ways to obtain the entanglement entropy Sa for a region A starting from the 
boundary and ending at r. First, we could consider the limit r% — > in S(ri,r), making the 
bulk region [r\,r] approach the boundary. As the twist field at r\ approaches the boundary, the 
correlation function diverges, because the presence of the boundary changes the regularisation 
necessary around the branch point. A way to evaluate the divergency is to use boundary 
conformal field theory, which applies in massive models when a local field is near to a boundary. 
It tells us that for small r\ there is a power law determined by the conformal dimension of 
T: fl(0|---T(ri)|0) B oc r7/ 2An [26]. We may then define T(0)|0) B as lim ri _ rl An T(n)\0) B - 
This appropriately regularised operator T(0) is simply proportional to the unitary operator 
performing a Z„ transformation, since its branch cut, through which 7L n transformations are 
performed, now extends through the whole space. But since |0) B is invariant under such a 
transformation, we find, with appropriate choice of proportionality constants, 

..boundary (r) = _ ^ ^ ^ ^ 

n-s-i an 

where we have changed the power of e in order to keep a scaling dimension of (essentially, this 
accounts for the change of regularisation necessary around the branch points). We used the fact 
that the entropy is real in order to change T — > T by complex conjugation. 

Second, we may take the limit r 2 — > oo in S(r,r 2 ). Then, the two-point function in (12. 4|) 
reduces to its disconnected part: 

B (0\f(r 2 )T(r)\0) B ~ B (0|t(oo)|0) BB (0|T(r)|0) B . (2.6) 

In the first factor, the twist field does not feel the presence of the boundary, hence this expec- 
tation value can be replaced by its expectation value in the model without boundary, (0|T|0). 
Dividing out this factor and using the appropriate branch-point regularisation, we find again 

Some of these considerations are made clearer by using crossing in order to implement the 
boundary as a state: 

B (0|f (r 2 )T(n)|0) B = (0|f (r 2 )T (r x )\B) . (2.7) 

The boundary state \B) is in the past at time in the Hilbert space of the model on the full line, 
and the twist fields are placed at imaginary times r\ and r 2 . More precisely, the boundary state 
is the n-fold tensor product of single-copy boundary states. The state (0| is the ground state of 
the n-copy model on the line, corresponding to asymptotic conditions at positive infinite times. 
No factor occurs in using crossing symmetry since the branch-point twist fields are spinless. The 
normalisation of the boundary state \B) is such that (0\B) = 1. Using crossing, we get 

5 boundary (r) = _ ± Zn£ 2A n ^ T ^ B y ^ 

n->i an 

It is this approach that we will be exploiting in this paper. 
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2.2 Large and short distance behaviour 



In the calculations above, we did not keep track of the normalisation constants occurring in 
reaching (|2.8f) . These correspond to an additive constant to the entanglement entropy, which 
is not universal. However, as mentioned in the introduction, the difference between constants 
occurring at large and short distances is universal. We now show that the expression (12, 8\\ gives 
the choice of large distance behaviour (|1.6p . When r is very large, we may use a decomposition 
of the identity, between the twist field and the boundary state, in energy eigenstates (see (|2.16|) ). 
The leading term comes from the ground state, giving 

^ oundary (r)~-lim-^-Z„ £ 2A "(0|T|0) (m^oo) (2.9) 
n-»l dn 

using (0\B) = 1. From [6], this is exactly — (c/6) log(em) + U/2 as defined in (|1.7p . hence this 
shows (|1.6p . 

We may give a CFT expression for the short-distance constant V(k). From the general 
formula for one-point functions in a boundary CFT with boundary state \B) [26j, we have 

(0\T(r)\B) ~ (T\B) CFT (2r)~ 2A " (2.10) 

as rm — > 0. Here, \T) is the normalised highest weight state corresponding to the primary field 
T in CFT. In the cases of finite perturbing boundary parameter h, the state \B) CFT is the UV 
limit h — > of the boundary state \B). In the case where h — > oo before rm — * 0, the state 
\B) CFT is the corresponding IR limit. Hence from (12.80 . we have 



5 boundary (r) _ £ log(2r / e ) _ H m ±( T \B) CFT + oil) {rm 0) (2.11) 
6 n-*i dn 

which gives 



V(k) = - lim ^-{T\B) CFT (2.12) 
n— >i dn 

where we recall that {0\B) CFT = 1, with (0|0) = 1 and (T\T) = 1. We see here that V{k) (with 
the definition (|1.5f) for k) takes only two values, one for k finite (UV) and one for k = — oo (IR), 
in agreement with (jl.8p . This is the main conclusion that we can derive from (|2.12p . as it is a 
non-trivial matter to evaluate this expression, and there may be subtleties associated to massive 
ground state degeneracies (see section [6]) . 

Note finally that we can write V(k) solely in terms of entanglement entropies: 

V(k) = lim (s h A °™ d ^(n) - S^^ir,- 1 ) - i^fa) + ^C^ 1 ) " fa*) , (^) 
where boundary and bulk entanglement entropies may be evaluated in different cut-off schemes. 



2.3 Integrable models and large-distance expansion 



For simplicity, let us consider a model of integrable QFT with a single particle spectrum and no 
bound states. The asymptotic states forming a basis of the Hilbert space are characterised by a 
number of particles k and by the rapidities 9j of these particles. As usual in QFT, two bases can 
be defined, representing particles coming in (m-states) and particles going out (out-states). In 
the n-copy model, we will denote the m-states by \6i, . . . ,0k)ai,...,a k with 0\ > ... > $k, where 
a,j are the copy labels; the vacuum will be denoted |0). The scattering matrix describes the 
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linear relations between the two bases. The two-particle scattering matrix of the n-copy model, 
that depends only on the rapidity difference 9 = 9\ — 62 by relativistic invariance, is given by 

S ab (9) = S(0) s " b , for a,b = l...n, (2.14) 

with S(9) being the scattering matrix of the original model and a, b the copy labels of the two 
particles. That is, the different copies of the model do not interact with each other. 

As mentioned above, the state \B) is just a tensor product of boundary states in the individual 
copies. In integrable models, these have an explicit expression as the famous boundary state 
introduced by Ghoshal and Zamolodchikov |22j . In the case where no boundary bound state 
can form, we have 

\B) = exp U-p^ jT * (y - 0) H~9)ZM j |0>. (2.15) 

The function R(9) is the boundary reflection matrix of the integrable QFT and Zj{9) are the 
Faddeev-Zamolodchikov operators, which provide a generalization of the creation-anhilation 
operators for integrable QFTs with non-trivial interactions |13|. I27j. Their main properties are 

Z ai (0i)...Z Ofc (0 fc )|O) = \6 U ---,0 k )a x ,...,a k ior9 1 >...>9 k 
Z a {8 l )Z h {8 2 ) = 5 a6 (0i - 9 2 )Z b (9 2 )Z a (9 1 ). 

The tensor-product form of the boundary state indicates that particles living in different copies 
of the theory do not interact through the presence of the boundary. 

One can now expand in (|2.8|) the boundary operator defined above to obtain a large-r 
expansion: 

00 

(0\T(r)\B) = (T)Y,fiVrm,K), (2.16) 



where 



1 " J- roo / 

(T)Mt, K ) = £ I]/ d9 r e~^ R ( 



VK 

T 



31 J2r"j'(=i 

FS hjlj2h --- jd \-9 l ,9 l ,...,-e i ,9e), (2.17) 



and 



Fp 1 -^, ...,e e ) = (oir^ZjM ■ ■ ■ z jt (e t )\o) (2.18) 

are the ^-particle form factors of the operator T in the n-copy model. For example: 

f (t,K) = 1, (2.19) 

h(t,n) = [°° d9R(^-9)F 2 Tlu (-9,9)e- tc ° shd J (2.20) 



4tt(T) 

and so on. This is a useful expansion, because the form factors of branch-point twist fields can 
be obtained exactly, in principle, by solving a set of consistency equations, first given in [6]. We 
will calculate these form factors in the next section for the Ising model. 
Employing (12. 8\\ and (II. 4|) . the entanglement entropy is given by 

c U °° 

SA(rm) = —- log(em) + — + st(2rm, k). ( 2 -21) 
d 2 L — ' 
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with 



Si(2rm, k) 



dfe(2rm,K) 
dn 



n=l 



(2.22) 



and where 




-4A. 



■71 



(2.23) 



is a universal constant, that relates the large-distance saturation to the short-distance logarith- 
mic behaviour in the bulk case. In obtaining the expression (|2.2ip for the entanglement entropy, 
we used the fact that T becomes the identity operator at n = 1, so that (T) n= x = 1 and all its 
form factors with one or more particles vanish. 

It is this form factor expansion that we will use in the sections that follow in order to evaluate 
the entanglement entropy in the boundary Ising model and in particular the constant V(k). 

2.4 Integrable boundaries in the Ising model 

The massive Ising model is characterised by the fact that the two-particle scattering matrix, 
as defined above, is S(9) = —1 (in the single-copy model): the particles are free Majorana 
Fermions. The central charge in (|2.2ip is 1/2 and the constant 



was obtained in j6l [28]. Let us now recall the types of integrable boundary conditions that 
have been found for the Ising model. A family that was studied in much detail in [22] is that 
corresponding to the presence of a magnetic field that couples to the Ising spin field (4>(t) in 
the action (|2.1|) ) on the boundary. To be precise, the spin field is the order parameter, hence 
we are looking at the scaling limit of the Ising spin chain in a transverse magnetic field whose 
magnitude is slightly below its critical value, and with a parallel magnetic field on the boundary. 
The corresponding boundary reflection matrix is given by 



which includes, for special values of the parameter k, the following physically different types of 
integrable boundary conditions, 

• Free boundary condition: n = 1, 

• Fixed boundary condition: k = — oo, 

• Magnetic boundary conditions (interpolating between the previous two): k = 1 — 
where h is a boundary magnetic field < h < oo. The free boundary condition would 
then correspond to h = whereas the fixed boundary condition is equivalent to having a 
infinitely large magnetic field fixed at the boundary. 

Boundary corrections to the expectation values of the energy and disorder field in the Ising 
theory were computed using this reflection matrix in [29|. 

In the cases where k > 0, the reflection matrix has a pole on the imaginary line on the 
physical sheet, < Im(#) < in /2. This implies that the boundary state expression (|2.15p is not 
correct. A modified expression exists [22] , but for simplicity, in the present paper we will not 
analyse this case. Hence, throughout we will consider k < (except in sub-section 13. 2p . Note 



U = -0.131984... 



(2.24) 
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that the case k = does not require modifications, since the residue of the i?-matrix vanishes at 
this point. At k = 0, the bound state becomes weakly bound, and propagates far into the bulk 
The cases k > —1, i.e. h < h c = 2y / m, are also somewhat special. In these cases, the 
R-matrix still has a pole on the imaginary 9 line, although not on the physical strip when k < 0. 
As noted in [22], the case k = —1 corresponds to a "critical" value of the magnetic field, h c , 
at which the reflection matrix happens to have a third order zero at 8 = 0. We will discuss 
the meaning of this critical point from the viewpoint of the boundary entanglement entropy in 
section 

Finally, we note also that the R matrix at k = is just equal to the negative of the fixed- 
boundary condition R matrix, k = — oo, and that 




(2.26) 



3 Form factor expansion for the entanglement entropy in the 
boundary Ising model 

3.1 Higher-particle form factors of Ising branch-point twist fields 

In order to evaluate the series (12.170 we need to address first the issue of computing higher- 
particle form factors of the twist fields in the Ising theory. In [6] the full set of form factor 
consistency equations was written but only the two-particle form factors were explicitly cal- 
culated. However, for the Ising theory, this will be the main piece of information needed, as 
higher-particle form factors can be obtained out of two-particle ones by using Wick's theorem. 

Since the branch-point twist field is invariant under the internal Z2 symmetry of the Ising 
model, characteristic of the Majorana Fermions, only even-particle form factors will be non-zero. 
Let us consider some of the consistency equations for the form factors of the twist field T (with 
even number of particles k) 

F^ ai - ak (e 1 + 27vi,...,e k ) = F^ a2 - ak(ai+1 \9 2 ,...,e k ,9 1 ) (3.2) 

where in we identify the particle types a + n = a. Using these relations repeatedly, it is possible 
to write all form factors in terms of form factors involving only one particle type (say 1). For 
1 < o,j < n and with the ordering 01 > . . . > a k , we have 

F T\ ai ...a k{e ^ _ A) = F p 1 - 1 (9 1 + 27ri( ai -l),...,9 k + 27Ti(a k - 1)), (3.3) 

and different orderings can be obtained using (13. ip . by which extra signs may appear. Using this 
as a definition for form factors with at least one particle of type different than 1, it is possible 
to check that equations (13. ip and (13. 2p are indeed satisfied, under the condition that (|3.ip holds 
for all particles being of type 1: 

-^fe"' 1 ' ' >0i>0i+l> ■••) = ~ F k ]{l l {- ■ ■ f0i+l>Qu ■ ■ •) (3-4) 
and under one additional condition, coming from n applications of (13.20 : 

Fp 1 -\9 1 + 2mn, ...,9 k ) = -F^ 1 - 1 (9 1 , ...,9 k ) (3.5) 
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where we used the fact that the number of particles is even. That is, the set of equations f)3. 1 j) 
and (13, 2h is consistent, and it is sufficient to solve (|3.4f) and (|3.5|) , 

Finally, there are two more conditions on form factors. We will write them in terms of form 
factors with only particles of type 1. These are the kinematic residue equations: 

- iRes- 6o=eo F^ 2 -\9 + ivr, 00,0! . . . A) = F fc T|1 - 1 (0 1 , ...,9 k ), (3.6) 



and 



- iRes^ifJ" 1 - 1 ^ + 2mn - nr,0 o ,0i . . . ,9 k ) = -F fe T|1 - 1 (0i, . . .,6 k ). (3.7) 

Since we are dealing the free Fermion case, it is natural to expect that the form factors of 
the twist field would admit closed expressions in terms of Pfaffians, as for the order and disorder 
fields of the Ising theory. This is indeed the case, and it is easy to show that 

F fc T|11 - 1 (0 1 ,...,0 fc ) = (T)Pf(K), (3.8) 

(recall that the Pfaffian has the property that Pf(iT) 2 = det(K)) where K is an anti-symmetric 
k x k matrix, k even, with entries 

Kij = r|1 f n 1 1 oJ — ■= K(9i - 9 3 ), (3.9) 

Tin,, f0\ /x 2nsinh(^)sinh(^) 
F T(0) = -i^-) and P{9) = liJ{*) ^ 

That is, these are the only non- vanishing form factors and are general solutions to (|3.4p -( f3~7"|) . 
Note that the Pfaffian expression is nothing else than the application of Wick's theorem on the 
operators Z\{9) involved in the form factors (12.180 (specialised to all particles being of type 1), 
a contraction of Z\{6\) with Z\{9<i) being K{9\ — 62). Hereafter the following properties of the 
function K{9) will often be used: 

K{9) = -K(-9), (3.11) 
K(9)\ n=1 = 0, (3.12) 
{K(9 + is)f = -K{9-is), 6,seR, (3.13) 

where "*" indicates complex conjugation. 

It is easy to prove that (|3.9|) solves (|3.4p - (|3.7p . Equation (|3.4p simply means that if we 
exchange two lines and the two corresponding (under the transpose) columns of a matrix, its 
Pfaffian gets a minus sign. Equation (|3.5p is a consequence of the property K{9 + 2iim) = —K{6) 
and (]3.1ip . and the fact that if we change the sign of a line and the corresponding column of 
a matrix, its Pfaffian also gets a minus sign. Finally, the kinematic residue equation (|3.6p can 
also be easily proved from the structure of the Pfaffian, seeing it in terms of Wick's theorem. 
Looking at the most singular term as 0o — * #o> using the poles of K(9) itself at 9 = ±iir with 
residues ±i, we obtain (|3.6p . 



3.2 Two- and four-particle boundary corrections to the entanglement entropy 
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3.2.1 Two-particle correction 



Employing (|3.8p for the two-particle form factor, it is not difficult to evaluate the first correction 
to the saturation value of the entanglement entropy (|2.21|) using (|2.22p and (|2.2U|) . We find 



si(t, k) 



k — cosh 9 J V cosh 2 9 



(3.14) 



since 



d_ 

dn 





y^\29) 


n=l d 


n 


iir tanh 8 


1 — 




dn 

n=l 


P{26)_ 


1 2 cosh# 



(3.15) 



The correction si(t, k) is finite for all values of t, including t = 0, as can be seen in figure [TJ At 
this point it is possible to evaluate the integral above explicitly: 



1 7r it v/1 — k (tt + 2arcsin(K)) 
ci(rc) := si(0, k) = h — ; — LL . 



and in particular 



. , 10-37T , . . 7T-2 

ci(— 1) = and ci(0) = — - — . 



(3.16) 
(3.17) 



In the case where k > 0, as was said in sub-section 12. 4} the boundary state expression (12.151) 
needs modifications due to the presence of a boundary bound state. Following |22j . the boundary 
state in the multi-copy model has the expansion 



\B) 



l + g^Z^+gl £ 2^(0)^(0) / ^i?( T -0)Z J H?)Z#) + 



d,e 



ITT 

2 



10) 



(3.18) 



where g c is proportional to the residue of the i?-matrix at the bound-state pole. Since the 
branch-point twist fields have zero one-particle form factors, the only possible modifications to 
the result for s\(t, k) above come from the quadratic term g 2 c Y^jk=i Zj(Q)Zk(.fy- The correction 
to si(t, k) would then be 



91^ 



d_ 

dn 



_d_ 

dn 



n=l 



il—l 

nY,K{2mj) 

j=0 



(3.19) 



n=l 



Using the methods of summation of [6], it is possible to show that this correction vanishes. 
Hence s\(t, k) is the correct leading large-distance correction for all k < 1. 



3.2.2 Four-particle correction 

Let us consider now the four-particle boundary correction. From (|2.17p and (|2.22p we find 

2 



S 2 (t,K) 



k e 



-t cosh 8 1. 



d_ 

dn 



E l rpT\iijj 
W) 4 



'1,171, — ^2, ^2) 



n=l 
(3.20) 
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where 



'1, Cli 



n-l 

«£ Fl ]111 \-0i,9i, -0 2 + 2m j, 8 2 + 27ui). (3.21) 



Here we have used (|3.ip and (|3.3|) . Employing (|3.8|) we find 

n-l 

n ^ ^1111 ^ _^ + 27r .^ ^ + 27rij) 
\ ' i=o 



(3.22) 



n-l 



n 



£ \K{29 1 )K(29 2 ) + K(9 12 - 2ma)K(9 12 + 2ma) + if (0 12 - 2-iria)K(-9 12 - 2-nia) 

a=0 

n—1 

K(29 1 )K(29 2 ) + n £ (#(012 - 2nia)K(9 12 + 2vrm) - if (0i 2 - 2iria)K(9 12 + 2vria)) . 



a=0 



Here and below we use 9; 



9j and 



9i + 0,-. It is simple to show that the 



n K(29\)K(29 2 ) term will give no contribution to the derivative at n = 1 (this is due to 
property (|3.12|) ) so that only the terms in the sum remain. These terms will give a contribution, 
since, employing (|3.13j) . they can actually be rewritten as 



n-l 



n 



£ [\K(9 12 - 2vrm)| 2 - \K{9 12 - 2maf 



a=0 



(3.23) 



Writing things in this way is useful as we can employ one of our main results in |6] , namely that 



d 

dn 



n-l 



n 



^2\K{9 -2ma)Y 



a=l 



(3.24) 



n=l 



Under integration in (|3.20p . the term with 9\ 2 can be changed into 9\ 2 , and the change of variable 
required inverts the sign of R(iit/2 — 9 2 ) using property (|2.26p . Hence both terms give the same 
contribution. Thus we find that the four-particle correction to the saturation value of the Ising 
entanglement entropy is given by 



s 2 (t, k) 



1 

32 



oo />oo 



-oo J — oo 
I roo 



'2 | 5(#i 2 )e- 2 * cosh ^r cosh ^ 



32 



-I d9R !E 



V 2 



-2t cosh t 



1 f°° 

— / d9 

32 



k + cosh 9 \ 1 — cosh ( 



(3.25) 



-It coshe 



k — cosh 9 J 1 + cosh ( 



Contrarily to the two-particle correction, we find that the four-particle contribution is divergent 
as t — ► 0. Technically, the reason for this is that the integrand of (|3.25p is a function that tends 
to the value — 1 as 9 —* oo when t = 0. Therefore the integral at t = is divergent. In order 
to find the precise behaviour of the correction as t approaches 0, one can rewrite the integral 
above as: 



dd 



k + cosh 9 \ 1 — cosh ( 



K — cosh 9 J 1 + cosh ( 



+ 1 



-2t cosh 6» 



^Ko(2t). 



(3.26) 
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The behaviour of the Bessel function as t goes to zero is well-known, 



K (2t) = -7 -log(t) + 0(t 2 logt), 



(3.27) 



where 7 = 0.577216... is the Euler-Mascheroni constant. Written in this form, the integral part 
is now a finite constant at t = 0, and we may define 



C2O) 



1 f°° 

— / de 

16 







K + cosh#\ 1 — cosh# 



cosh# 



+ 1 



1 



(1 + k) 5 



with in particular 



2k - 3k 2 



<*(-!) 



1 + 



23 
120' 



1 + cosh 9 
k (2k — 1) (it + 2 arcsin At) 



vT 



(3.28) 



(3.29) 



The constant C2(k) is an increasing function of the magnetic field (i.e. a decreasing function of 
k). Therefore we have 

1 , . 7 



s 2 (t,n) 



16 l0g(t) + T6 +C2(K) + 0W 



(3.30) 



(the order 0(t) comes from the next correction to the integral part of (|3.26p ). Note that in 
view of the short-distance behaviour of the entanglement entropy (|1,6|) . we expect that the 
coefficients of the logarithmic divergencies at small rm will add up to the finite number 1/12 
when all corrections are considered. This will be proven in section [5j Also, the constant c^^k), 
like c\{k) above, is a part of the constant V(k) in (|1.6p : again, in principle one should add up 
all such constants, for all corrections, in order to obtain V(k). This will be discussed also in 
section [5j 

An interesting mathematical phenomenon can be observed: although the integral in (|3.26p 
has the same value for k = —00 as for k = for any t > 0, we have C2(— 00) = 3/8 different from 
c 2(0) = 1/8. The explanation is that the limit t — > of the integral in (13.260 as a function of k 
is not uniform. For all values of t > we have s^ft, —00) = S2(t, 0), and there is a maximum for 
k E (—00, 1) at a unique value k = Kq. But as t becomes smaller, the position of this maximum 
shifts towards more negative values, until it reaches —00 at t = 0. There, if we take away the 
constant (as function of k) term ^ log(i) in order to make the limit finite, the value of the 
maximum itself reaches C2(— 00). It is also possible to observe in the integral in (|3.28p that 
the symmetry between n = — 00 and k = is broken. Indeed, if k is very negative, the term 
in parenthesis can be approximated by 1 except for values of 9 where k + cosh 9 ~ 0. These 
are very large values of 9, but they are not damped by any other factor, hence the mistake in 
approximating by 1 is non- negligible for any n. 

This means that formally, the expansion (]3.30p is valid only for k > —00. For the case 
k = — 00, that is, the fixed boundary condition, we have to consider the other order of the 
limits: first k — > —00, then t — > 0. By the symmetry between k = —00 and k = 0, we define 



c 2 (-oo) =: c 2 (0) = - 



(3.31) 



so that (|3.30p still holds in the case of a fixed boundary condition, k = —00. 

In fact, it is instructive to obtain a more general small-t expansion, where we take simulta- 
neously k — * —00. Let us consider t — > with — nt = a fixed. We may use the change of variable 
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s = cosh# — 1 and write S2(t, k) as 



16 y \k-l-s) V (s + 2)V2 + s + l) e 16 J Q ds { k _ 1 _ s ) s + 1 ■ 

(3.32) 

The first integral as a function of n has a uniform limit as t — * on k £ [— oo, 0), so that we can 
directly take k = — oo and f = 0; this gives (2 — log2)/16. The second integral does not have a 
uniform limit, but it can be evaluated explicitly: 

Uae- 2t - (a + t)T(0, 2t)e 2t - 8a(a + t)e 2(a+t) r(0, 2(o + *)) 



16(o + t) 

where T(z, u) is the incomplete Gamma function, v z ~ l e~'"dv. The small-t limit can then 
easily be taken: 

s 2 (t, -a/t) = 4 log(t) + ^ + 4(a) + 0(t). (3.33) 
lb lb 

where 

4(a) = jj-iae 2a r(0,2a). (3-34) 

It is easy to see that c\{a) interpolates between lim^-oo C2(re) at a = to C2(— oo) at a = oo. 

We expect this phenomenon of non-commutativity of the limits to be generic: it should occur 
at all orders, except, as we have seen, for the very first two-particle correction. Its meaning will 
be explained in sectionEJ We also refer the reader to section[6]for an analysis of the large-distance 
corrections found here. 



3.3 Higher particle boundary corrections to the entanglement entropy 

As one would expect, the form factor expressions which are obtained from (|3.8p become more 
and more involved as the number of particles is increased. In particular, it is easy to see that 
the 2£-particle form factor of T is made out of the sum of 

(21 - IV 

1 x 3 x 5 x ... x {21 - 1) = 2 i 1{£ _ J 1)V (3.35) 

terms. This number grows faster than exponentially with I. However, as we have already seen for 
the 4-particle corrections, not all of these terms contribute to the derivative at n = 1 and those 
that contribute give in many cases the same contribution due to the fact that they are identical 
when integrated over. For example, the two contributions in (13.231) are actually equivalent when 
integrated over 9\ and 02- The combination of these two factors, that is, terms whose derivative 
vanishes at n = 1 and terms that can be grouped together, allows us to reduce very dramatically 
the amount of non-vanishing contributions to the derivative that we obtain from higher particle 
form factors. 

The 2£-particle form factors contributing to the entanglement entropy (see (|2.22p and (|2.17p ) 
occur in a sum of the form 

n 

F^^ 1 '"^^ e (—9i,6i, —02,02, ■ ■ ■ , —0i, 0e) 

h,-,je=o 

n-l 

= " E FS l -\-9 u e l , (-9 2 y\9i\. . . , (-Oiy'-^e?- 1 ), (3.36) 

3l,—dt-l=0 
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where 

e b a = 9 a + 2mb. (3.37) 

Here, we used (|3.3|) and (|3.1|) . along with the fact that every copy number occurs in pairs, so 
that no sign remains. Prom Wick's theorem, one of the terms contributing to this sum will be 
the following contraction: 



? T|l...l, 

2i 



(-ex,0i, (-0 2 y\ 0i\ (-0 s ) j2 , ■ ■ ■ , &£h(-oi-iY t -\#£i(-o t y'-\# l i - x ), (3.38) 



which corresponds to 



n-l 



(3.39) 



We will call such a term "fully connected" , which denotes any contraction where the sums over 
j's cannot be factorised into a product of sums. 

Thanks to (|3.12j) . any fully connected term vanishes as n — > 1. Below we will show that 
the derivative with respect to n is not zero as n — ► 1, but converges to a distribution in the 
rapidities, generalising the main result (|3.24p of [BJ. This means that the product of two or 
more fully connected term has a vanishing derivative as n —* 1, so that the only terms that will 
contribute to (|3.36|) are the fully connected Wick contractions. We will further show below that 
all fully connected terms can be brought to the form (13.380 . 



3.3.1 Explicit evaluation of fully connected terms 



The sum (13.390 can be obtained in a systematic way by exploiting a result which was first 
obtained in appendix 3 of [6j. The result derived there was a special case of the sum 

" — \ ?'sinh ( y ~t x } 

(3.40) 

s cosn — <•( >si i — 

a=0 



2 cosh I cosh I 



which can be computed in exactly the same manner. The sum (|3.39p can be evaluated by simply 
using (13,400 recursively. When doing so one realizes the need to distinguish two special cases, 
depending on whether t is even or odd in (|3.39|) . The final expressions are, 



n-l 



E K ({- xl r)K(4^) . . . Kixitr^K^r: 

3l,—,32i-l=0 



(-l)^sinh( ^f Xl 



21 

11 2 cosh 
i=i 



21 -\ 
1-3 



2i: 



2t 



if (E * + (2j - l)ni) + K(Y^ Xi - (2j - 1)™) 



i=i 



i=i 



(3.41) 



and 



n-l 



(-l) m 2cosh( 



l^i=l ■ L 



3i,— ,he=0 



2t 

i 



v 21+1 I 

wE^+E 

' t=i j=i 



21 
i-j 



2i+l 

II 2 cosh f 
i=i 

2t+l 2i+l 

K{^Xi + 2mj) +K{Y j x 1 - 2nij) 



i=l 



i=l 



(3.42) 



A proof by induction of these identities is provided in appendix [A"l 
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3.3.2 Analytic continuation in n and computation of the derivative at n = 1 of fully 
connected terms 



Let us start by considering the analytic continuation and evaluating the derivative at n = 1 of 
(|3.4ip . This summation formula already provides an analytic continuation in re of the multiple 
sum over j's for any fix x's. Naturally, such an analytic continuation is not unique. Additionally, 
this formula itself gives many more analytic continuations when taken under integration over 
the x variables. Indeed, for integer re, all poles of the functions K involved are exactly cancelled 

by the zeros of sinh( ^ — i ), hence the integration contours can be moved away from the real 
a>axis without any change to the answer. But for non-integer n, poles and zeros generically are 
at different points, except for the poles of K at values of its argument ±i-7r. Hence, different 
integration contours give different analytic continuations of the integrals. 

Below we suggest that the formula (|3.4ip is the correct one, but that the contours need to 
be taken differently. Let us consider the rapidity variables 6>j, with xi = (and x^i = 02t,i)- 

For n large enough, all the poles of K are beyond the integration contour of the variable 



for all j's in (|3.4ip . However, as n decreases, poles cross the 6 integration contour. At the 
points where these poles cross, there is a zero provided by the hyperbolic sine function, so 
that no discontinuity occurs. Yet as function of n, the result is not smooth, since there is a 
discontinuity in the derivative with respect to n. This is because as re approaches an integer 
value where a pole crosses the integration contour, the right-hand side of (|3.4ip is not uniformly 
convergent as function of 8, and develops an infinitely thin peak of finite height at 6 = 0. There 
is no natural way of modifying (|3.4ip in order to avoid this phenomenon. Indeed, it is related 
to the fact that for integer values of n, the unambiguous value of the sum in (|3.4ip at 8 = 
varies with n up to n = £ + 1, but from n = £ + 1 up to infinity it takes the same values. Hence, 
any analytic continuation will have to reproduce this unnatural behaviour. We note that this 
phenomenon is a generalisation of that observed in [B] in the two-particle approximation of the 
bulk two-point function (and above in the four-particle boundary case): there the value of the 
sum at 6 = at n = 2,3,4,... is constant and non-zero, while the value at n = 1 is zero. 
The resulting non-uniform convergence was at the basis of the calculation of the derivative with 
respect to n at re = 1. 

In order to recover a smooth function of re up to n = 1, we then need to move the 6 
integration contour towards values where the argument of K is just 6 ± iir in (|3.4ip . for all j's. 
It is convenient to still avoid poles of the hyperbolic cosine factors in the denominator. Hence, 
we consider shifting all 9i by ^^-ni, so that we get the following summation formula, valid 
under integration: 

£ K({-9 l2 )^)K{6^ K@%?) L 2*sinh0 



h,— ,he-i=o 



11 

K{26 + qm) \[ shift. 



J=H=± V 7 n2cosh(%i-^7ri 

i=l v 

Here, shifty . )g . j- 1 ^ is an operator acting on all other functions in the integrand (that is, 
those not appearing here), indicating that their arguments must be modified as written. More 
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precisely, this indicates a shift of contour, and if the other functions in the integrand have poles 
that are crossed by this shift, then the residues must be taken. It is important, of course, that 
the integrals over rapidities stay convergent at all stages of the shifts. 

The strongest evidence for the validity of this formula is provided by the fact that in the 
bulk case, the scaling dimension of the twist fields can be reproduced by re-summing the form 
factor expansion, for any real n > 1. In appendix [O we provide very convincing numerics for 
this scaling dimension. 

We can now evaluate the derivative with respect to n at n = 1 of this. For generic 9, the 
function on the right-hand side of (j3.44H is zero at n = 1, but at 8 = 0, it is non-zero. Hence, 
we need to properly take the limit n — ► 1 of the derivative of the right-hand side of (|3.44|) as a 
distribution. We have 

-£[**«(* ± ^U = T I=™, (3.45, 

but this has a pole at 8 = 0, coming from the kinematic pole of K at ±(2n — l)iri. It can be 
resolved by noticing that for 8 — > and n — > 1, we have 



— [smh8K(28±iri) 
an 



7T 



2 6>T (n- 1)tu' 



(3.46) 



Hence, as a distribution, 



- — [smh8K(28±7ri)] n=1 



7T 



cosh ( 



2 sinh(0 =F «0+) 



and this gives 



(3.47) 



dn 



n-l 
31,— J2i-l=0 



iri cosh 9 



j=l q=± 



21 -\ 
L-3 



n=l 

21 

csch(6l - qi0+) ]J shift 6 
i=i 



21 

J! 2 cosh 

8=1 



j-1 
I 21 



Til 



(3.48) 



Finally, we may simplify this formula by shifting back the contours towards their initial positions. 
Doing so, the integrand will be zero, and the only contributions will be poles taken on the way. 
The final result is 



d 

dn 



E 

ii>— ,ht-i=0 



n=l 



j = l k=l q=± 



1/2 (k = j) 
1 (k < j) 



21-1 



2i 

n shift, 



21 

n 

i=i 



11 2cosh(^i+ (? ^7ri 



(3.49) 



-5(8). 



The case k = j in this formula can be simplified using 



27r 2 5(8) £ 



2£-l 



2£- 2 
i-l 



{-l)V5{8) 



2t 

n 2 cosh 



2£ 

]7J 2 cosh ,, 



(3.50) 
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so that we obtain 



d_ 

dn 



n-l 

E 

jl vj2£-l=0 



(3.51) 



J n=i 



-1) £ 27T 2 J(0) 



2£ 



1 



2r 



2 cosh , 
Lj=i 



33 + 1 



e j-i 

EEE 

j=l h=l q=± 



2£ 



J 



L \ . 2« 

(-1V f] shift, ft ^i- k 

1=1 ^ 



21 

n 2 cosh 



21 



Note that the first term inside the square brackets on the right-hand side is what is obtained 
by directly using (|3.4ip . without contour shifts; it is a direct generalisation of the four-particle 
case. The other terms are corrections, characteristic of higher-particle contributions only. 
A similar analysis can be made for the odd case, (|3.42p . The result is 



d_ 

dn 



n-l 

m-9i2T)K{& 

ji>—,he=0 



23 j 



K 027+l,l) 



I 



(3.52) 



j=l k=l q=± 



21 



J 



n=l 
21+1 

(-iy q n shift 6 

i=l 



. j-k+1/2 . 
i+1 21+1 m 



21+1 

]1 2 cosh 

8=1 



+ q 



j-fc+l/2 
2£+l 



TTl 



Note that in this case, we get "pure correction terms", as directly taking the derivative with 
respect to n at n = 1 of (|3.42[) . without contour shifts, gives zero. 

The final results as written in (|3.5ip and (13.521) hold only if the other functions of the 
integrand do not have poles on the region covered by the shifts; then the shift operators are just 
shift of arguments. The contributions of poles may be evaluated in a similar way, but below 
we will avoid these complications. Also, these formulae are valid for £ = 1,2,3,... only, that is, 
excluding £ = 0. 



3.3.3 Putting everything together 

The final step in order to evaluate se(t, k) is to work out the "multiplicity" of (|3.39p for a fixed 
particle number. This can be easily done by exploiting the Wick contraction picture introduced 
in (|3.38|) . Let us pick the first rapidity —9\. There are 2(£ — 1) possible contractions that could 
be performed as a contraction with 9\ is not allowed (that would produce a non-fully connected 
term). Let us assume that —9\ is connected to 9\ for some fixed i,j. Then, if we now pick the 
next rapidity, that is 9±, it can be connected to almost any term, except for —9\ (that would 
again be non-fully connected) and the two that are already connected. That gives us 2(£ — 2) 
possibilities. We carry on this argument by connecting 9\ to some 6>ij,, then 9\, to some 9 3 
etc., until no rapidities are left. Hence we find that the total number of fully connected terms 
is (£ — 1)!2^ _1 . All these terms are identical when integrated in all rapidities because they can 
all be brought to the form (|3.39[) by a series of two types of operation. First, we may change 
the sign of a rapidity, getting a minus sign from the R matrix thanks to property (I2.26p . then 
change the order of the two Z operators associated to this rapidity, which cancels this minus 
sign. Second, we may move pairs of Z operators associated to a given rapidity without getting 
any sign. 
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In order to be able to use formulae (|3.5ip and (|3.52p , we must make sure that no poles occur 
in the regions covered by the shifts, for any £. This imposes that the other functions in the 
integrand should not have poles for lm(0j) G (— 7r/2, 7r/2). Since factors R(iir/2 — Qj) occur in 
the integrand, we must take 

« < 0. (3.53) 



With this restriction, we obtain 
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both formulae hold for £ = 1,2,3,.... All integrals involved are absolutely convergent, and it 
is easy to see that this holds throughout all shifts required. Notice that for £ = 1 in (|3.54|) we 
recover the result (I3.25P as it should be. Setting £ = in (|3.55|) does not give (I3.14|) . since £ = 
is out of the range of applicability of these formulae; the two-particle case s\(t, k) is a special 
case. 

Hence, we have for the entanglement entropy in the boundary case 



S 



boundary 



(rm) 



~ log(em) + — + ^2 s e( 2rm i K ) 



(3.56) 



with (13~T41) . (133411 and (|3~55j) . 

Formulae (|3.54p and (|3.55p are quite lengthy, but can be written in a more symmetric way, 
more appropriate for numerical calculations. We present such alternative expressions in appendix 

m 



4 Form factor expansion for the bulk entanglement entropy 

We have studied the entropy of the Ising model in the presence of a boundary. However the form 
factor approach employed so far has first been used for bulk theories in [6] where the entropy 
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was obtained in the two-particle approximation. We should mention here that the entanglement 
entropy of free theories was obtained by a different method in [30, 31J. There it was shown how 
it is connected to Painleve transcendents. The analysis of the previous sections can be easily 
adapted now to find closed expressions for all entropy contributions in the bulk case. Employing 
a form factor expansion for the two-point functions of the twist fields in the bulk theory we can 
write 
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Since we want to compute the derivative at n = 1 of the above, we need to evaluate 
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where we used ()3.13p . 

The only contribution to the derivative (|4.2[) will again come from the fully connected terms 
for similar reasons as in the boundary case. In fact, it is not difficult to convince oneself that the 
non- vanishing contributions coming from the 2-Oparticle form factor will be completely analogous 
to the non-vanishing contributions coming from the 4-Oparticle form factor in the boundary case, 
the only difference being the presence of the reflection matrices in the latter case and the amount 
of terms that contribute. Let us consider a very particular fully connected term, corresponding 
to the Wick contractions 
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The contribution of this term to (|4.2f) is 
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By changing the sign of rapidities with odd index, which does not change the result under 
integration in (|4.ip . we obtain exactly (13.51H . 

As for the boundary case, it is possible to argue that all fully connected terms are identical 
to the one above when integrated in all rapidities. The number of such terms can easily be 
evaluated by the following argument: Let us consider the first form factor in the product f)3.35|) 
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and pick one of the rapidity variables on which it depends, say 9j for some fixed It can 
be connected to any other rapidity in the same form factor, say (P., , and therefore there are 
21 — 1 possibilities. Now we can look at the second form factor, picking 9Z J . There are 21 — 2 

possible connections, since it cannot be connected to QT 3 as this would produce a factorisable 

_jti 

term. Suppose it is connected to 9-,r ■ We now come back to the first form factor, looking 
i" 

at 6K„ . It can be connected to any other available rapidities, so there are 21 — 3. Continuing, 
we find that there are (2£ — 1)! fully connected terms. Under the integral all these terms are 
equivalent, because exchanging two rapidities in both factors simultaneously does not bring out 
any sign. Therefore, 

, CO 
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so that 
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Again, see appendix [B] for an alternative expression of e<(rm). In particular, the 1 = 1 contri- 
bution is given by 
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which is one of the main results obtained in [BJ. 
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5 Exact UV behaviour of the entanglement entropy in the bound- 
ary Ising model 

5.1 Exact logarithmic behaviour 



Let us start by extracting the exact small-rm logarithmic behaviour of the bulk entanglement 
entropy from the full form factor expansion (14. 7|) . For this purpose, after integrating Q21 using 
the delta-function, it is convenient to change variables to the set Xi = Oi^+i for i = l,...,2£-2, 
and #2^-1- Then, we have 
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and in particular, 
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Hence, the hyperbolic cosine factors in the denominators, involving all Xj's as well as the last 
two sums above, do not depend on 621-1, so that it is the large |#2£-i| behaviour of the integrand 
that determines the singularity at small rm. In the exponential, we have at large | #2^-1 1 an d fi x 
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Hence, the exponential factor will display a jump of finite width from 1 to around |#2^-i| ~ 
— log(rm) + const, when rm is small, so that the leading small-rm behaviour can be obtained 
by omitting the exponential factor and replacing J d02e-i by — 21og(rm). As a result, we find 



ei(rm) ~ ^log(rm) 
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where we used the fact that the integrals factorise into even and odd-indexed x variables. The 
function Ji(a) can be evaluated exactly: 
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The expression for the number gg can be simplified to 
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(5.9) 

The sum over i of the first term in (|5.8|) is readily seen to be 1/4, whereas the rest was verified 
to be consistent with —1/12 by summing 500 terms. Hence, we find 



(5.10) 
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in agreement with the known logarithmic behaviour from CFT. In appendix we provide a 
numerical analysis of the scaling dimension A n characterising the short-distance behaviour of 
the full two-point function of twist fields from our form factor expansion. Agreement with the 
known CFT dimension provides an extremely non-trivial check of the validity of this form factor 
expansion. 

5.2 Exact expression for V(k) 

From the results in the previous two sections, it is possible to derive an exact expression for 
V{k) appearing in (jl .61) . 

First, comparing (|4.6p to (|3.54[) . we see that if we set the reflection matrices to 1, we have 
2s2i (2rm, k) = ei(2rm). The meaning of this is that the replacement R — ► 1 in Y2i=i s 2i{2rm, k) 
precisely provides the leading small-distance behaviour of the entanglement entropy in the 
boundary case, which is related to the bulk case by S^° undary (rm) ~ ^S^ ulk (2rm) (rm — > 0) (up 
to a finite term). The term si(2rm, k) in the boundary case only contributes a finite term as 
rm — ► 0, and this holds as well as for all other terms S2£+i(2rm, k) for £ = 1, 2, 3, . . .. 

In order to understand this, consider first the expression for S21 in (|3.54p . The leading small- i 
behaviour of the integrals in the square brackets comes from the region of large rapidities which 
is not damped by the hyperbolic cosine factors in the denominator. Integrating 621 using the 
delta-function, we can change variables to the set Xi = ^j+i for i = 1, . . . , 21 — 2, and Bn—\. 
Then, we have 

21-2 

Oi = Y {-1?-% + {-i) 1+i e 2i -i (5.11) 

j=i 

and in particular, 

21-1 l-i 21-2 t-l 21-1 l-i 

Y 0i = X^i-i + 0m-i. Y 61 = Z^^i-i, Yl 9i = J2 X2 r ( 5 - 12 ) 

i-l j=l j=l j=l i=2 j=2 

Hence, the hyperbolic cosine factors in the denominators, involving Xj's and the last two sums 
above, do not depend on 621-1, so that it is the large |#2£-i| behaviour of the integrand that 
determines the singularity at small t. On the other hand, 02t-i is involved in the argument 
of every i?-matrix. Since R(in/2 — 9) ~ ±i as Re(#) — > ^00 for k finite, and the opposite for 
k = —00, one can see that the product of R matrices goes to 1 as #2^-1 ~> ice for any k. This 
shows that the replacement R — * 1 gives the leading behaviour as rm — ► 0. 

Second, let us consider a similar change of variable in the odd case (|3.55p . Zj = 9i,i+i for 
% = 1, . . . , 21 — 1, and 621- Then, 

2&-X t 21 t-i 

Y 0i = Y x 2j-i - 821, Y 6i = ^2 X2 i + e<2t ' ( 5 - 13 ^ 

i=l j=l i=2 j=l 

so that all variables are involved in the hyperbolic cosine factors in the denominators. Hence, 
in this case the limit t — > can be taken, and all integrals are still convergent. 

Then, subtracting this leading behaviour, and using the known small-rm behaviour of the en- 
tanglement entropy in the bulk case, (jl.7p . we find that the small-rm behaviour in the boundary 
case is given by (ll.6p with 
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with c\{k) given by f)3. 16|) and where 
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both formulae for £ = 1, 2, 3, . . .. These generalise C2(k) given in (I3.28|) (recall that 9 = Y^j Oj) 
The derivation of (|5.14j) is as follows: 
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Although we have provided arguments suggesting that the integrals defining q(k) are convergent, 
in order for (15.141) to be a correct representation of V(k), the infinite sum over £ should give 
a finite result. This is much more subtle, as form factor expansions are expected to provide 
convergent series expansion for finite distances, but not necessarily at zero distance. Exact 
evaluations of the first few coefficients q(k) for some k below, and extrapolation to higher £, give 
strong indications that the series is indeed convergent; hence (|5.14p is a correct representation. 

5.3 Exact and approximate re-summations 

All integrals defining the coefficients can be evaluated exactly. For instance, a change 

of variable Xj = e d i makes all integrands rational functions, and multiple integrals can bring 
logarithmic terms, which can all be integrated. We have been able to evaluate exactly all 
coefficients q(0): 
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extrapolating from the exact values at I = 1, 2, 3 and 4 obtained with the help of Mathematica. 
Along with ci(0) from (|3. 16j) . these re-sum to: 
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^From this, we directly obtain V(— oo), since C2i(— oo) = C2i(0) and C2£+i(— oo) = — C2^+i(0). 
This gives 
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We have also computed the values of q(— 1) for I up to 8 and obtained the following exact 
expressions, 
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These values already give a rather good approximation to V{— 1): 
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0.321966... 
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which is compatible with the expected value of log\/2 = 0.346574.... Therefore, as expected, in 
the UV limit the entropy only depends on whether the magnetic field h is finite (free boundary 
conditions) or infinite (fixed boundary conditions). An interpretation of these results for V(n) 
is provided in section El 



6 Discussion 

We recall that our main results are 1) the full form factor expansion for the entanglement entropy 
in the boundary Ising model (|3.56p for k < 0, 2) the full form factor expansion in the bulk case 
(I4.7|) . and 3) an expression for the short-distance constant V(k) (|5.14p leading to its exact value 
(jl.8p . We now discuss our results for the boundary entanglement entropy in the large-distance 
and short-distance limits. 
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6.1 Infrared behaviour of the boundary entanglement entropy 



The form factor expansion (|3,56p converges rapidly for large distances. Hence we expect the 
2- and 4-particle contributions evaluated in section (|3.2p to provide a good approximation of 
the IR behaviour of the entanglement entropy for sufficiently large distances. Let us start by 
studying the behaviour of the correction (I3.14p for t = 2rm large. In this case it is possible to 
compute the leading contribution of the correction by expanding the integrand (except for the 
exponential term) in the variable cosh — 1 around the value zero. We find, 

Sl(t,K / -1,1) 

and 

si(M) 

si(t,-l) 

Recall that these hold for all k < 1. 



1 PjT K + 1 e * f r/o, 

= — \ f + 0(e"*t~ 5/2 ) 




Figure 1: The function si(2rm, k) for several values of k. 
A similar computation for S2(t, k) yields: 

«<«■ ^ = -k$0$ £m + °™- 5/2 >. < 6 - 3 > 

and 

S2 U _i) = _ /I e " 2 * + 0( e - 2 '(2t)- 9 / 2 ), (6.4) 
A ' ; 1024V 2 (2i) 3 V2i 
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Figure 2: The function s\(2rm, k) + S2(2rm) for < r < 2. The value of rm for which the 
entropy is maximal for a given value of k increases as k gets closer to 1. As in figures [2] (a) 
and (b), the maximum of the entropy seems most pronounced for n = 0.7 and rm around 0.1 
and tends to shift and eventually disappear as k — > 1. However, since the maximum occurs 
generally for rather small values of rm, its position may be affected significantly by higher order 
form factor corrections. Moreover, for k > 0, additional contributions to S2 may come from the 
boundary bound state. Yet, it is certain that there is a maximum for k > —1 as in this case the 
function tends to zero from above as rm — > do. 

which is clearly sub-leading comparing to the s\(t, k) expansion. Note that there will be correc- 
tions to this sub-leading behaviour for k > due to the boundary bound state. 

Let us examine more carefully (16. ID and (16. 2p . The leading contribution to the entropy is 
negative, as long as k < —1: the asymptotic value is approached from below. For k > — 1, 
however, a change of sign occurs: the asymptotic value is then approached from above, as can 
be seen in figure [TJ This means that for values of k > — 1 the "saturation" value of the entropy 
in the infrared limit, U/2, is actually not the maximum value the entropy reaches. For k > — 1 
there exists some finite value of rm for which the entropy has a maximum. From looking at the 
curves of si(t, k) + S2{t, k), this value seems to be moving towards higher rm as k is decreased 
towards —1 - see figures [2] (a) and [2] (b). However, for —1 < k < the maximum is at very 
small values of rm, and for k > 0, S2 may receive extra boundary bound state contributions. 
Hence, this behaviour of the maximum cannot be conclusive. 

The fact that the entanglement entropy is not monotonic for k > — 1 is not in contradiction 
with its fundamental properties. In particular, the fact that the entanglement entropy is an 
increasing function of rm in the bulk case follows from the "strong subadditivity theorem" 
and translation invariance, as proven in |32} I33j. Since translation invariance is broken in the 
boundary theory, the entanglement entropy in this case is not necessarily a monotonic function 
of rm. 

As recalled in sub-section 12. 4[ for k > — 1, the R-matrix has a pole on the imaginary line, 
although it is only on the physical strip for k > (where a bound state is present). The value 
k = — 1 (corresponding to a magnetic field h c = 2^/m) is a "critical" value where the i?-matrix 
has a third-order zero at = 0. What we see here is that it is this critical value that plays an 
important role for seeing a maximum from the large-distance behaviour. Note that a similar 
sort of large-distance behaviour was observed for the one-point functions of the energy and spin 
operators in the boundary Ising model [29]. The explanation suggested was that while at short 
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distances one sees the free boundary condition for any finite h, at large distances one observes 
the fixed boundary condition for h > h c , and the free one for h < h c . This gives rise to a 
cross-over behaviour in the h-mr plane for h < h c , where we observe two separated regions, 
corresponding to increasing and decreasing entanglement entropy with distance. 

6.2 Ultraviolet behaviour of the boundary entanglement entropy 

6.2.1 Interpretation of V(k) 

It is natural to interpret V(k) as a "boundary entanglement": the contribution of the boundary 
to the entanglement between the region A and the rest. Naturally, for fixed boundary condition, 
there should be no contribution at all, since the boundary does not experience quantum fluctu- 
ations. Our result (|1.8jl shows that we have chosen the correct large-distance normalisation to 
have V{— oo) = 0. On the other hand, for free boundary conditions, the boundary fluctuates and 
should participate to the entanglement. This is in agreement with V(k > — oo) = log\/2 > 0. 

In fact, we may connect V(k) to the boundary entropy s, a quantity that essentially counts 
the number of degrees of freedom pertaining to a boundary. This quantity is simply given 
by s = log<7 where g is the boundary degeneracy introduced by Affleck and Ludwig [25]. In 
particular for a bulk CFT, they showed that g = (0\B) where |0) is the bulk CFT ground state, 
and \B) is a normalised boundary state in the bulk CFT Hilbert space (in particular, s < 0). 
It is considering the boundary entropy that Friedan and Konechny [3l] were able to provide a 
proof of the "g-theorem" : that the (/-function decreases in the RG flow from UV to IR. 

In order to derive the relation between V(k) and s, we first provide a direct connection 
between V(k) and entanglement entropies of the Ising model at criticality: 

V(k) = 5^ oundary (r) critical - ^ ulk (2r) critical + log y/2. (6.5) 

The result is independent of r, since both entanglement entropies have the same logarithmic 
r-dependence. Most importantly, both entanglement entropies must be evaluated in the same 
cut-off scheme; for instance, both should be evaluated on the same lattice, with the same lattice 
spacing. The entanglement entropies should be evaluated at the UV fixed point if k is finite, 
and at the IR fixed point if k = —00. 

In quantum field theory, it is not easy to implement the same lattice spacing in the bulk 
and boundary situations, since the lattice spacing enters into non-universal constants. Using 
massive QFT, it is possible to solve this problem. Consider the entanglement entropy SU( r i> ?"2) 
of sub-section 12.11 but with a slightly different normalisation specified below; we will denote it 
Sa(ti, T2). We may uniquely fix the cutoff, for instance by requiring the conformal normalisation 
as above, iSaC^Ij^) ~ f l°g(( r 2 — r i)/ £ ) + as r 2 — * T\. This, then, is just the bulk critical 
entanglement entropy above: 5 , ^ ulk (r) critica i = | log(r/e). But using the same object, hence with 
the same lattice spacing, we can also define the critical entanglement entropy in the boundary 
case, 5 , ^ oundary (r) CI . itica i, following the arguments already outlined in sub-section 12.11 Note first 
that for r2 3> ri ^> 0, the entanglement entropy saturates to some constant — |log(me) + U 
thanks to the presence of the mass, and this saturation is a sum of the contributions of the 
two boundary points at n and ri- These contributions are equal, so that one boundary point 
contributes — | log(me) + U/2. Now for r2 3> r\ and r\m finite, we get the entanglement entropy 
in the boundary case (i.e. with the region ending on the boundary being connected) and with 
a mass, but with an extra contribution of the boundary point r^ (at infinity). Hence, we may 
define 5^° undary (r) = lim^^oo S(r, ^2) + ! log(me) — U /2. This quantity has large-rm asymptotic 
given by —5 log(me) + U /2. We then only need to take the critical bulk limit rm — > 0. 
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Figure 3: How a vacuum of the form • • • <g) |0, +) <8> |0, — ) (g) |0, +) (g> • • • becomes a two-particle 
state after going through a branch-point twist field. 

At this point it would seem natural to identify U with U, the bulk saturation constant, 
and then have S^° undary (r) = S^ oundary (r). However, this does not hold. The main idea comes 
from the fact that in the bulk, the ground state is degenerate, with spins being asymptotically 
up or down, |0, +) and |0, — ). If we allow fluctuations amongst the two degenerate ground 
states, as would be obtained from a large-volume limit of a periodic space, the bulk constant 
U should be identified with the two contributions from the boundary points, and an additional 
pure entropy contribution coming from this fluctuation, which should be log 2 (for instance, an 
extra contribution of log 2 due to zero-modes was present in results of |35j). Hence we must 
write 

U = U + log 2. (6.6) 

In the boundary case with any non-zero magnetic field, however, there is a definite ground state 
that is chosen at infinity, and the derivation above indeed involves only the boundary point 
contribution U/2. Hence the large distance limit of 5^° undary (r) is — |log(me) + U/2 — log \[2 

,i , r ifTTTi riboundary / \ ^boundary / \ i fK mi n j riboundary , n 

so that from (|1.6|) . b A \r) = b A (r) — logy 2. Ihen we nnd b A (^critical = 
| log(r/e) + V(k) — log y/2, which shows (|6.5j) . 

Technically, allowing fluctuations amongst the two degenerate ground states corresponds to 
choosing |0) = ((|0, +) + |0, — ))/v2)® n in the n-copy model. Our branch-point twist field form 
factors as constructed in [6], and more generally here, assume the use of this ground state, since 
they are non-zero even in cases where there are sheets with odd numbers of particles. Recall that 
in the ordered regime, particles correspond to domain walls, separating spin-up and spin-down 
regions. Then, the only way to have out-states with, for instance, two particles on different 
sheets after the twist field, is to have an in-state with sheets on different ground states, see 
figure[3j A completely up ground state |0, +)® n , for instance, wouldn't admit such form factors. 
Since all configurations of an even number of particles distributed amongst the n sheets are 
allowed, we must use the completely symmetric ground state. 

The constant U is the large-distance limit of the bulk entanglement entropy evaluated from 
(0,+|® n T(ri)T(r 2 )|0,+)® ri , whereas U is that obtained from (0|T(n)T(r 2 )|0) with the sym- 
metric ground state. The difference can be computed explicitly. First note that matrix elements 

(0, e' x \ ® ■ ■ ■ ® (0, e' n \f (n)T(r 2 )|0, ei) <g> • • • ® |0, e n ) 

have zero large-distance limit unless = ej = e' for all i,j (that is, all signs are the same), since 
otherwise domain walls will have to propagate between the twist fields. Hence we immediately 
find 

(0|r(ri)T(r 2 )|0) ~ 2 1 ""(0, +|® n t(r 1 )T(r 2 )|0, +)® n (6.7) 
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at large distances, so that, taking derivatives with respect to n, we have 

- (A (0 |f( ri )T(r 2 )|0>) _ ~log2- (J-(0, +|^(ri)T(r 2 )|0,+>®») _ (6.8) 
which gives (16. 6|) , 

Note that although the arguments above hold for finite, and perhaps only large enough, 
magnetic field, the constant V(k) is the same for any k finite (see sub-section 12. 2p . hence the 
result (|6.5p is valid in general!. 

Finally, Calabrese and Cardy [36] proposed a formula for the boundary entropy, which 
amounts, from (|6.5p . to the statement 

V(k) = a + log y/2, (6.9) 

and which has been tested numerically in |37[ [38] This is in agreement with our results (|1.8|) , 
since s = in the free boundary case, and s = — log\/2 in the fixed boundary case [26]. 

From these arguments, in general the shift log v2 should be replaced by — log C where C 2 is 
the fraction of the ground state degeneracy broken by the boundary condition for large enough 
h. This implies (11.90 . Such an "extra" contribution to s = log g was also found in [39] in 
calculations using Thermodynamic Bethe Ansatz techniques and was accounted for by means of 
similar arguments (it is important to note that the Thermodynamic Bethe Ansatz used there is 
fundamentally different from our approach based on form factors). 

In [39] the flow of g between critical points was studied for several families of minimal Toda 
field theories with ground state degeneracy k and with a boundary completely breaking it. The 
factor C = 1/ \fk was termed "symmetry factor" . Interestingly, it was shown that at the infrared 
point of these models, one gets g = C. Hence, for massive models with spontaneously broken 
order-parameter symmetry, and with an order-parameter boundary perturbation, one should 
find V(k) = s — sir > 0, where sir = logC is the infrared value of s. For the same order- 
parameter perturbation both on the bulk and boundary, one should simply find V{k) = s < 0. 
All these considerations should not depend on integrability. 

Note that a more detailed study of V(k) may be useful in understanding if the boundary 
entropy s is always bounded from below. 

6.2.2 Conformal bulk with non-conformal boundary 

Taking the short-distance limit mr — ► while increasing the magnetic field with a fixed product 
2nmr = a, we obtain a conformal bulk theory with a non-conformal boundary, with a = h 2 r the 
dimensionless parameter relating the distance r to the boundary and the boundary magnetic field 
h. The entanglement entropy then interpolates between the UV free-boundary point r — > 0, and 
the IR fixed-boundary point r — > oo. We observed explicitly this interpolation in the four-particle 
result in section 13.21 through the monotonic function c\ (a) (|3.34p . It is a complicated matter to 
generalise this to higher particle contributions, but it would be interesting to understand if the 
function stays monotonic, as this gives a natural ^-function flow. 

7 Conclusion 

The main technical results of this paper are exact infinite-series formulae for the bi-partite 
entanglement entropy of the Ising model with and without boundaries. We have used for this 

*We would like to thank here P. Calabrese and P. E. Dorey, for suggesting to us the possibility of a relationship 
between the degeneracy of the ground state and the log \/2 term in (|6.5|l . 
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a relationship between entanglement entropy and the derivative with respect to n at n = 1 of 
correlation functions of branch-point twist fields in the model composed out of n non-interacting 
copies of the Ising model. In order to obtain our formulae, it has been necessary to tackle several 
non-trivial intermediate problems: 

• finding closed expressions for all non-vanishing form factors of branch point twist fields in 
the n-copy theory, 

• identifying the correct analytic continuation in n of the contributions of these form factors 
to correlation functions and evaluate their derivatives with respect to n, and 

• checking both the form factor formulae and their analytic continuation for consistency. 

The first of these problems was easy to solve, as we are dealing with a free theory for which 
Wick's theorem applies. Hence, we have been able to show that all form factors of the twist 
field admit expressions in terms of a Pfaffian. Obtaining the right analytic continuation of every 
contribution to the form factor expansion of the twist-field boundary one-point function (or 
bulk two-point function) is a highly non-trivial problem. The main complication arises from the 
fact that the pole- and zero-structure of the form factors (in particular, the way some poles and 
zeroes cancel each other) changes substantially as soon as n is allowed to take non-integer values. 
As a consequence, the phenomenon observed in [6] of non-uniform convergence of form factors 
as n — ► 1 is generalised to a non-uniform convergence as n — > £' for all positive integers £' < I 
for the 2£-particle form factor. The problem was solved following the principle that the analytic 
continuation is obtained from the analytic function that describes form factor contributions at 
values of n large enough. This amounts to evaluating first the contribution we would obtain 
analytically continuing in n around some integer and then adding the residues of all the extra 
poles that are crossed by the integration contours when bringing n from infinity. This is however 
not the only analytic continuation that is possible and therefore it becomes quite crucial to find 
ways of checking it for consistency. In this paper we have been able to do this in a very precise 
manner by finding an explicit formula for the leading logarithmic behaviour both of the two- 
point function of the twist field in the bulk and of its derivative at n = 1 (that is, the bulk 
entanglement entropy). By a combination of analytical and numerical computations we have 
been able to extract both UV behaviours with extreme accuracy and to show that they agree 
with what is expected from CFT arguments. 

The results just described have put us in the position to analyse another quantity of interest 
in this context, that is the contribution to the free energy that can be attributed exclusively to 
the presence of the boundary. This is essentially the boundary entropy, the natural logarithm of 
the boundary degeneracy or ^-factor originally introduced by Affleck and Ludwig in |25j . In our 
analysis we have computed the universal quantity V(k) which is closely related to g (where n is 
related to the boundary magnetic field in the Ising model). We have defined V{k) as a certain 
rm independent contribution to the boundary entanglement entropy in the UV limit. We have 
found an exact formula for V(k) and evaluated it exactly at k = 0, — oo to V(0) = log\/2 and 
V{— oo) = 0. We have also gathered strong numerical evidence that V(k) is in fact constant and 
equal to V(0) for any finite values of k. These two values of V(«) would correspond to the two 
conformal invariant boundary conditions that are known for the Ising model: the free and fixed 
boundary conditions. The fact that V(k) is larger for free boundary conditions (finite magnetic 
field) and that V(0) — V(— oo) = log\/2 are properties which also hold for log g. However, it is 
known from Cardy and Lewellen's work that gf ree = 1 and gfi xe d = 1/V% hence 



V(k) — log g = log\/2. 
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In this paper we have identified the difference between these two quantities as an IR contribution 
to the entanglement entropy coming from the ground state degeneracy in the periodic case that 
is broken in the boundary case. We have proposed a generalisation of this to more general 
models with relevant boundary perturbations. 

We have also found interesting results for the IR behaviour of the entanglement entropy: it 
is known that the bulk entropy saturates for large distances (in the Ising model, the saturation 
value is given by the constant U given in (|2.24p ). In particular, the fact that the entropy is 
an increasing function of rm follows from the "strong subadditivity theorem" and translation 
invariance, as proven in |32} 133]. Since translation invariance is broken in the boundary theory, 
the entanglement entropy in this case is not necessarily a monotonic function of rm. Indeed, we 
find that, for a range of values of k, it has a maximum for some value of rm before reaching its 
asymptotic value U/2. This range of values of k starts precisely at k = —1 which corresponds 
to a value of the magnetic field for which, in a sense, the boundary becomes "critical" . 

There are many open problems related to the present work and in general, to the computation 
of the entanglement entropy in integrable QFT. In the case of the bulk Ising model, it is known 
that the entanglement entropy can be described via Painleve transcendents [301 EI] ■ It would be 
interesting to check the consistency of this representation with our full form factor expansion. 
In the boundary Ising case, it is known that the one-point function of the order parameter has 
a Fredholm determinant representation for any magnetic field, from which differential equations 
can be derived |29j . It would be interesting to see if similar formulae hold for the branch-point 
twist field. In the general QFT case, the most obvious problem is perhaps to extend the present 
analysis to theories other than the Ising model. We believe that this should be possible to some 
extent but it is very unlikely that re-summations can be done analytically for interacting models. 
Yet an independent check of (|1.9p in the more general situation would be useful. It would also be 
interesting to apply the form factor approach employed here and in [8] to the computation 
of the entanglement entropy of multiply connected regions, both for bulk and boundary theories 
and also to extend the analysis of the present paper to the finite temperature situation. 

Acknowledgments: We are very grateful to P. Calabrese, J. L. Cardy and P. E. Dorey for 
discussions and interest in this work. We also would like to thank the Galileo Galilei Institute for 
Theoretical Physics in Florence for hospitality during the program "Low-dimensional Quantum 
Field Theories and Applications" where key discussions occurred. B.D. thanks the Mathematics 
Group at City University London and O.C.A. thanks the Mathematics Department of Durham 
University, where parts of this work were completed. O.C.A. acknowledges financial support 
from a City University London "Pump Priming" fellowship. 

A Re- summation of fully connected terms 

In this appendix we would like to provide a proof by induction of the equalities (I3.4ip - (l3,42p . 
Let us then assume that (|3.41|) is true for some value of t and try to obtain (|3.42[) from it. If 
(|3.4ip holds, then it will also hold when the variable xn is shifted as X2t — ► X2t — 2nip. We will 
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do this and introduce at the same time one extra sum in the variable p and one factor K(y p ), 

E m- xl r) K {4-») . . . k^-^k^^kw = tp^Sl 

h,-,ht-i,P=o n 2 cosh 

8=1 

n-l i 

< 

where, as before x := X)i=i x i- We can now em pl°y (|3-40p in order to carry out the sum in p in 
the second line, 



n-l l , op _ -, x 

E E ( 7_ / J t^"" - ivr ) + ^ (*~ J '~ P + K (y p ), (A.i) 



^ z cosh ( ~ 

E ^ (* ±j - p T ^)^(y p ) = 2sinh x CQ 2 sh y (^0* + 2/ ± 2 ^') +K(x + y± 2vri(j - 1))) . (A.2) 

p=0 2 2 

Substituting these sums into fO) we obtain 

E «-< 4 i m4 • • • Ki^tr^Ki^r^w) = 1 J z 2 c i osh( — ) 

h,-,h£-i,P=0 2 cosh | f| 2 cosh ^ 

i=l 

x E ( 2 /j 1 ) + y) + + v) + K &~ j + y) + K( x - j+1 + y )] . (a.3) 

i=i ^ J ' 
The sum in j can be split as, 

e-l 



E ( 2 /_ ) ) + f ) + + f )] + E ( e 2 - i - 1 ) + v) + ^ + v)] 

j = l \ J / j=Q \ J / 

* )K(x + y)+J2[ i U 3 ) + V) + + y)] - ( A - 4 ) 



where we have used the identities 

_ / 21 - 1 \ / 2£ \ , / 2£-l\ ( 2t-\ \ ( 21 \ , , 

This completes our proof, since by calling y = X2£+i our expression (IA.3j) is nothing but (|3.42jl . 
In an entirely analogous way it is possible to obtain (|3.4ip starting with (|3.42p . 



B Alternative formulae 

In this section, we re- write formulae (|3.54|) . (|3.55|) . (|4.6|) . (|5.15|) and (|5.16|) . These are obtained 
using 5(9) = f 2^e t>ie and shifting all contours in 9 back to the real line. The results involve 
less and more symmetric multiple integrals, which is useful for numerical computations or exact 
evaluation of the integrals with the help of a symbolic mathematics computer application (like 
Maple or Mathematica) . 
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we have 
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C Scaling dimension from the form factor expansion of the bulk 
two-point function 

The scaling dimension of branch-point twist fields was evaluated from their two-particle form 
factors, along with the known form factors of the trace of the stress-energy tensor, for arbitrary 
n > 1 in [6 J using the A-sum rule [ID]. This gave an exact result (in agreement with the CFT 
result), because in free models, the stress-energy tensor does not have form factors with more 
than two particles. The scaling dimension can also be recovered from the two-point function 
of twist fields themselves, by analysing its short-distance behaviour. This is not as convenient, 
however, since it involves re-summing an infinite number of terms, which contain more and 
more integrals (all form-factor contributions), and usually it is superfluous, since the A-sum 
rule already gives the scaling dimension. But in our case, it is cornerstone check, because the 
correct analytic continuations in n of higher-particle contributions to the two-point function of 
twist fields necessitated a non-trivial choice of integration contours, a problem which did not 
occur for the two-particle contribution. 

The most convenient way to evaluate the scaling dimension from the form factor expansion 
of the two-point function in free models is to take the logarithm, so that we are looking for the 
coefficient of the small distance logarithmic divergency. Using the fact that all form factors can 
be evaluated using Wick's theorem, the logarithm of the two-point function is just the sum of all 
connected contributions. The analytic continuation in n of connected form factor contributions 
was obtained in (|3.44p . Using the same symmetry arguments as in section H] in order to count 
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the number of connected contributions of a given particle number, we find 



( (0\f(r)T(0)\0) \ ^ 1 
8 1 (T) 2 I ^o, 



2( 



21 

n 



_ j —rm cosh 0j 

2tt 



(C.l) 



where ui{9\, . . . , 9<it) i s n times the right-hand side of (|3.44|) . 

As in sub-section 15. Q it is possible to obtain the small-rm logarithmic divergency by inte- 
grating out the 02t variable, after the change of variables to the set Xj = Ojj+i, j = 1, • • • , 21 — 1 
and 02£. Separating into even and odd-indexed x variables and using Fourier transform, we can 
evaluate all remaining integrals in a similar fashion to what was done in section [57TI In fact, it 
is convenient to shift back all contours in (|3.44p to unshifted contours, picking up residues. The 
contribution to un{9\, . . . , 9 2 i) of the sum of all unshifted contours is simply 



n{— l) l 2i sinh 9 



E 



21-1 
i-j 



Fi(9) 



(C.2) 



where 9 = X)ff=i @j an< ^ 

F t (d) = K{29 + (2j - l)m) + K(29 - (2j - l)wr) 
For the coefficient of — log rm this gives a contribution 
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The sum of the residues taken when unshifting the contours gives the following corrections to 
this contribution: 
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where [•] means the integer part, and the function J^(a) was evaluated in (15. 7p . Here, the sum 
over q is restricted as follows: 



E' 

q=± 



for k = [j/n], the term with q = +1 is present if and only if n[j/n] < j — 1. 



The sum of (|C.4p and (|CJ .6|) should be compared with 4A n = (n— l/n)/12. 

Taking about 30 terms in the sums over I in both (IC.4j) and (|C.6|) and adding up both 
contributions gives the following numerical values, for n between 1 and 3 (see fig. U|): 
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Figure 4: The exponent in the short-distance behaviour of the bulk two-point function, as 
function of the number of sheets n analytically continued to n > 1. This should be twice the 
scaling dimension, and the curve shows the expected value from CFT. The data points come 
from a re-summation of about 30 terms in the form factor expansion for the logarithm of the 
two-point function. For n > 2, an improvement using Euler's formula is needed in order to make 
the series convergent. 

The corrections (IC.6|) are exactly zero when n is an integer: in these cases, no contour shift 
is necessary from the beginning. However, the numerical evaluation of (|C4p shows that this 
contribution is insufficient to reproduce the correct dimension formula for non-integer n, showing 
the necessity of the contour shifts. It is also interesting to note that the sum over I in (|C,4p is 
in fact a divergent alternating sum for n > 2. One can however make it convergent by using 
Euler's formula for improving convergence of alternating sums: 

OO OO I , _ v 

5> = ^2-V = fc-i W (C7) 

e=i i=\ k=i ^ ' 

The rationale behind this is that the sum is in fact convergent for any rm > 0, where Euler's 
formula can be used, which then gives a finite limit as rm — * 0. 
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